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A tagged medium-energy neutron beam has been used in a precise measurement of the absolute 
differential cross section for np back-scattering. The results resolve significant discrepancies within 
the np database concerning the angular dependence in this regime. The experiment has determined 
the absolute normalization with ±1.5% uncertainty, suitable to verify constraints of supposedly 
comparable precision that arise from the rest of the database in partial wave analyses. The analysis 
procedures, especially those associated with the evaluation of systematic errors in the experiment, 
are described in detail so that systematic uncertainties may be included in a reasonable way in 
subsequent partial wave analysis fits incorporating the present results. 

PACS numbers: 13.75.Cs, 21.30.-x, 25.10.+S 



I. INTRODUCTION 



Theoretical treatments and applications of the 
nucleon-nucleon (NN) force at low and intermediate 
energies have progressed considerably in sophistication 
through the past decade. Partial wave analyses and po- 
tential model fits to the NN scattering database have 
incorporated explicit allowance for breaking of isospin 
(/) symmetry, e.g., by removing constraints that previ- 
ously required equal I — 1 phase shifts for the pp and 
np systems, and have been used to constrain the pion- 
nucleon-nucleon coupling constant . Effective field the- 
ory approaches have become competitive with more 
traditional meson-exchange models of the interaction, in 
terms of the quality of fit provided to the database and 
the number of adjustable parameters employed, while 
holding out the promise of providing internally consis- 
tent two- and three-nucleon forces from the same theory. 
Striking success has been achieved in ab initio calcula- 
tions of the structure of light nuclei 0] by combining phe- 
nomenological three-nucleon forces with NN interactions 



•Present address:Dept. of Physics, Texas A & M University, Col- 
lege Station, TX, USA 

tPresent address: Dept. of Radiology and Radiological Sciences, 
Vanderbilt University, Nashville, TN, USA 

-^Present address: Dept. of Physics, University of Zagreb, Zagreb, 
Croatia 

^Present address: Dept. of Physics, George Washington University, 
Washington, DC, USA 



taken without modification from fits to the NN scatter- 
ing database. An important aspect in these advances has 
been the approach toward consensus on which measure- 
ments should be included in an NN database to which 
conventional optimization techniques can be sensibly 
applied. The rejection of specific, allegedly flawed, ex- 
periments from the database has not been without con- 
troversy. In the present paper, we report detailed results 
from a new np scattering experiment addressing one of 
the most prominent of these controversies. 

Discrepancies among different experiments have led 
to a drastic pruning of cross section measurements for 
intermediate-energy np scattering. For example, the 
SAID partial wave analysis (PWA) of the np database 
01 rejects more than 40% of all measured cross sections 
in the range 100-300 MeV in neutron laboratory kinetic 
energy. The rejected fraction is even larger in the Ni- 
jmegen PWA especially so for scattering at center- 

of-mass angles beyond 90°. The rejected data include 
nearly all of the most recent experiments, carried out by 
groups at Uppsala and Freiburg . The problems 
are illustrated in Fig. ^ by the comparison of data from 
these two groups with earlier Los Alamos measurements 
|lOl | that dominate the medium-energy back-angle cross 
section data retained in the database. Clear differences 
among these data sets are seen in the shape of the angu- 
lar distribution. Other differences, reflecting the general 
experimental difficulty in determining the absolute scale 
for neutron-induced cross sections, are masked in the fig- 
ure by renormalization factors that have been applied 
in the partial wave analyses. Removal of the Uppsala 
and Freiburg data, which exhibit fairly similar angular 
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dependences, begs the question of whether the crite- 
rion used to reject them 0, IMS 01 subtly bias the 
PWA resuhs toward agreement with older measurements 
that might have had their own unrecognized systematic 
errors. 

The np back-angle cross section discrepancies have 
been highlighted in debates concerning the value and ex- 
traction methods for the charged ttNN coupling constant 
(in the notation of pseudovector formulations of the 
interaction, or equivalently g'^^/Air in pseudoscalar for- 
mulations) 0,01 . np scattering PWA's appear to deter- 
mine this basic parameter of the NN force well: e.g.^ the 
Nijmegen analysis 5j yields = 0.0748 ± 0.0003 (equiv- 
alent to g"^ — 13.54 ± 0.05), and the authors claim 
that the constraints are imposed by the entire database, 
with no particularly enhanced sensitivity to any specific 
observable In contrast, Ericson et al. have extracted 
a significantly higher coupling constant, consistent with 
older "textbook" values {g'^/An w 14.4), by applying con- 
troversial pole extrapolation techniques to the Uppsala 
back-angle np scattering cross sections alone 0|. While 
much debate has centered on the rigor of the pole extrap- 
olation method 0, Q, 0| , it is clear that the discrep- 
ancy in coupling constant values arises in large part [ig 
from the cross section discrepancies between the Uppsala 
measurements and the "accepted" database. An exper- 
imental resolution of these discrepancies is highly desir- 
able, especially if a new experiment can also pin down the 
absolute cross section scale. Bugg and Machleidt have 
pointed out 0] that the largest uncertainty in their de- 
termination of is associated with the normalization 
of np differential cross sections, which are often allowed 
to float from the claimed normalization in individual ex- 
periments by 10% or more in PWA's. In contrast, the 
Nijmegen group claims that, despite sizable normal- 
ization uncertainties in existing elastic scattering data, 
precise total cross section measurements fix the np ab- 
solute cross section scale to ±0.5% accuracy. This claim 
could also be checked by a new experiment that provides 
good experimental precision on absolute differential cross 
sections. 

In the present paper we report detailed results from 
such a new experiment, designed to resolve these np back- 
angle cross section discrepancies. The experiment em- 
ployed techniques completely independent of those used 
in other medium-energy measurements, in order to pro- 
vide tight control over systematic errors. A tagged neu- 
tron beam [l3| centered around 194 MeV kinetic energy 
bombarded carefully matched, large-volume CH2 and C 
targets, which permitted accurate subtraction of back- 
grounds from quasifree scattering and other sources. The 
bombarding energy range was chosen to match approx- 
imately that used in earlier high-precision np scattering 
polarization measurements from the Indiana University 
Cyclotron Facility (lUCF) [11 E ll^. Recoil protons 
from np scattering were identified in a detector array 
of sufficient angular coverage to measure the differential 
cross section at all cm. angles beyond 90° simultane- 



ously. 

The tagging allows accurate determination of the ab- 
solute scattering probability for the analyzed subset of 
all neutrons incident on a secondary target, but it also 
offers a host of other, less obvious, advantages important 
to a precise experiment: (1) accurate relative normaliza- 
tion of data taken with CH2 vs. C targets; (2) event- 
by-event determination of neutron energy, impact point 
and incidence angle on the secondary target, with the 
latter measurement being especially important for cross 
section measurements very near 180° cm. scattering an- 
gle; (3) three-dimensional location of background sources 
displaced from the secondary target 0| ; (4) precise mea- 
surement of the detector acceptance for np scattering 
events; (5) methods to tag np scattering event subsam- 
ples that should yield identical cross section results but 
different sensitivity to various sources of systematic error. 
The tagging was thus essential to the entire approach of 
the experiment; no extra work was required to extract 
absolute cross sections, and thereby to provide an im- 
portant calibration standard for medium-energy neutron- 
induced reactions. 

The basic results of this experiment have recently been 
reported brieffy [23|. In the present paper we provide 
more detail on the comparison of results to PWA's, on the 
data analysis procedures and on the evaluation and char- 
acterization of systematic uncertainties. Such details are 
important for resolving the sort of discrepancies that have 
plagued the np database. Partial wave analyses should, 
in principle, incorporate experimental systematic as well 
as statistical errors in optimizing fits to data from a wide 
variety of experiments. In order to do so, they must have 
access to clear delineations of which errors affect only the 
overall normalization, which have angle-dependence and, 
in the latter case, what the degree of correlation is among 
errors at different angles. Overall systematic uncertain- 
ties in the absolute cross sections reported here average 
±1.6%, with a slight angle-dependence detailed herein. 
Statistical uncertainties in the measurements are in the 
range ±(1 — 3)% in each of 15 angle bins. 



II. DETAILS OF THE EXPERIMENT 

The experiment was carried out in the lUCF Cooler 
ring l22l| , with apparatus (see Fig. ^ installed in a ring 
section where the primary stored proton beam was bent 
by 6°. A primary electron-cooled unpolarized proton 
beam of 202.5 MeV kinetic energy and typical circulat- 
ing current of 1-2 mA was stored in the ring. Neutrons 
of 185-197 MeV were produced via the charge- exchange 
reaction p+d n±2p when the proton beam passed 
through a windowless internal deuterium gas jet target 
(GJT) of typical thickness « 3 x 10"^^ atoms/cm^. The 
ultra-thin target permitted detection of the two associ- 
ated low-energy recoil protons from the production re- 
action in double-sided silicon strip detectors (DSSD's) 
comprising the "tagger". Measurements of energy, ar- 
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FIG. 1: Comparison of previous np scattering differential cross section measurements (a) from Uppsala 73 and Los Alamos 
and (b) from PSI |P| and Los Alamos (T^ near 200 MeV. The Los Alamos data in each case are represented by closed 
squares and the other data by open circles. The experimental results are compared to the Nijmegen PWA93 ^ partial wave 
analysis solution evaluated at appropriate energies. The Los Alamos data have been renormalized by factors of 1.092 in (a) 
and 1.078 in (b) to bring them into agreement with the PWA. The relative cross sections reported in have been similarly 
normalized here, while the reported absolute cross section scale for the Uppsala data has been retained. 



rival time and two-dimensional position for both recoil 
protons in the tagger, when combined with the precise 
knowledge of cooled primary proton beam direction and 
energy, allowed four-momentum determination for each 
tagged neutron on an event-by-event basis. During the 
measurement periods, the stored proton beam was oper- 
ated in "coasting" mode, with rf bunching turned off to 
minimize the ratio of accidental to real two-proton co- 
incidences in the tagger. The proton beam energy was 
then maintained by velocity matching (induced naturally 
by mutual electromagnetic interactions) to the coUinear 
electron beam in the beam cooling section of the rmg. 

Details of the layout, design and performance of the 
tagger detectors and of the forward detector array used 
to view np scatteringevents from the secondary target 
are provided in Ref. 17]. Here, we summarize the salient 
features briefly. The tagger included an array of four 
6.4x6.4 cm^ DSSD's with 480 fj,m readout pitch in two 
orthogonal {x',y') directions, each followed by a silicon 
pad ("backing") detector (BD) of the same area. The 
DSSD's were positioned about 10 cm away from the cen- 
ter of the gas jet production target. Each DSSD had 
128 x' and 128 y' readout channels. Readout was ac- 
complished with front-end application-specific integrated 
circuits (ASIC's) that provided both timing and energy 
information [T^ . The timing signals provided to exter- 
nal electronics consisted of the logical OR over groups of 
32 adjacent channels of leading-edge discriminator sig- 
nals based on fast shaped and amplified analog signals 
generated in the ASIC's. The timing signals available 
from 4 a;' x Ay' logical pixels for each DSSD permitted 



operation of the tagger in a self-triggering mode, where 
the time-consuming digitization of slow pulse height sig- 
nals from all 1024 DSSD channels could be initiated by 
logic based solely on the tagger hit pattern, as recon- 
structed from the fast timing signals. This self-triggering 
was critical to the determination of precise absolute cross 
sections, because it allowed acquisition of data to count 
directly the flux of tagged neutrons that did not interact 
in the secondary target or in any of the forward detectors. 

Only recoil protons that stopped either in the DSSD's 
{Ep < 7 MeV) or BD's {Ep < 11 MeV) were considered in 
the data analysis, because for these the tagger provided 
a measurement of total kinetic energy with good reso- 
lution. By combining these energy measurements with 
position measurements for both recoil protons, we were 
able to determine the energy and angle of each tagged 
neutron, within their broad distributions, with respec- 
tive resolutions of as ~ 60 keV and (Jangle ~ 2 mrad. As 
part of this determination, we reconstructed the longi- 
tudinal origin (zvertex) of each produced neutron within 
the extended GJT density profile with a resolution of « 
2 mm, by comparing neutron momentum magnitudes in- 
ferred by applying energy conservation (independent of 
vertex position) vs. vector momentum conservation (de- 
pendent on vertex position) to the tagger information for 
the two recoil protons. Similar resolution was obtained 
for the transverse coordinates at which each tagged neu- 
tron impinged upon a secondary scattering target (TGT 
in Fig.O positioned 1.1 m downstream of the GJT. 

Two solid secondary targets were used during the pro- 
duction running: a 20 x 20 x 2.5 cm"^ slab of polyethy- 
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FIG. 2: Top view of the np scattering experiment setup. 



lene (CH2) containing 1.99 x 10^^ hydrogen atoms/cm^ 
and a graphite target of known density machined to have 
identical transverse dimensions and the same number of 
carbon atoms per unit area. Each target thickness was 
determined to ±0.4% by weighing. Data were collected 
in 18-hour cycles, comprising 6 hours of running with 
the CH2 target, followed by 6 hours with C and 6 more 
hours with CH2. The frequent interchange of the tar- 
gets facilitated accurate background subtractions. Both 
targets intercepted neutrons over an approximate pro- 
duction angle range of 14° ± 5°, and cuts were generally 
placed on the tagger information during data analysis to 
confine attention to tagged neutrons that would hit the 
secondary target. Such tagged neutrons were produced 
at a typical rate of ~ 200 s~^, leading to typical free np 
back-scattering (angle-integrated) rates ~ 1 s^^ from the 
CH2 target. 

Protons emerging from the secondary target were de- 
tected in a forward array of plastic scintillators for trig- 
gering and energy information, and a set of (three-plane) 
multi-wire proportional chambers (MWPCs) for track- 
ing, as indicated in Fig. |21 The plastic scintillators in- 
cluded large upstream veto (LUV) and small upstream 
veto (SUV) counters to reject charged particles produced 
upstream of the secondary target. The AE scintillator 
was separated from the secondary target by a MWPC to 
permit easy discrimination against np scattering events 
initiated in that scintillator. The rear hodoscope com- 
prised 20 plastic scintillator bars jl^ of sufficient thick- 
ness (20 cm) to stop 200 MeV protons and give 15-20% 
detection efficiency for 100-200 MeV neutrons. All for- 
ward detectors were rectangular in transverse profile, 
with the rear MWPC and hodoscope spanning a con- 
siderably larger vertical than horizontal acceptance. The 
entire forward array provided essentially 100% (> 50%) 
geometric acceptance for np scattering events initiated at 
the CH2 target for angles 6'c.,„. > 130° (Ocm. > 90°). For 
cm. angles forward of 90° the large size and significant 



neutron detection efficiency of the hodoscope provided a 
small efficiency for detecting forward-scattered neutrons 
in coincidence with larger-angle protons that fired at least 
the first two MWPCs. 



The tagger and forward detector array were designed 
to facilitate a kinematically complete double-scattering 
experiment with a first target giving negligible energy 
loss. With the same apparatus, a similar measurement 
of pp scattering was possible simultaneously. For this 
purpose one could use the tagger to detect a single large- 
angle recoil deuteron instead of two recoil protons, in or- 
der to tag a secondary proton beam via pd elastic scat- 
tering in the GJT. By requiring a coincidence between 
a single hit in the tagger and a signal from the small 
upstream veto scintillator (SUV in Fig. 12), we could de- 
fine a secondary proton beam of very similar transverse 
dimensions to the tagged neutron beam. Another scintil- 
lator (Veto2) placed just in front of the rear hodoscope 
allowed us to distinguish, at trigger level, between pro- 
tons from pd elastic scattering that traversed the forward 
array without further nuclear interactions and protons 
that scattered out of this secondary beam in material 
following SUV. In the present paper, we discuss only the 
former group, as their yield provides an accurate relative 
normalization of runs taken with the CH2 vs. C targets. 



The triggered events of interest for the present analysis 
were recorded in four mutually exclusive event streams, 
three for tagged neutron candidates (consistent with two 
distinct tagger hits and no accompanying signals from 
LUV or SUV) and one for tagged proton candidates (con- 
sistent with a single tagger hit in prompt coincidence 
with both LUV and SUV). The trigger logic defined these 
event streams as follows: (1) tagged neutrons with no 
rear hodoscope coincidence, providing a prescaled (by a 
factor of 20) sample for neutron fiux monitoring; (2) np 
scattering candidates for which a tagged neutron was in 
coincidence with signals from both the AE scintillator 
and the rear hodoscope; (3) tagged neutrons in coinci- 
dence with the hodoscope but not with the AE scintilla- 
tor, a sample used for evaluating the neutron detection 
efficiency of the hodoscope 0; (4) tagged protons in 
coincidence with both the AE and Veto2 scintillators, 
providing a prescaled (by a factor of 10) sample includ- 
ing pd elastic scattering events from the GJT, used to 
cross-normalize C and CH2 secondary target runs. The 
total flux of tagged neutrons intercepting the secondary 
target was determined from a sum over event streams 
(l)+(2)+(3), while comparative analyses of the three 
streams facilitated crosschecks to calibrate the system 
and aid understanding of potential systematic errors. 
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III. DATA ANALYSIS 

A. Cuts and conditions on tagged neutron beam 
properties 

The general philosophy of the data analysis was to de- 
fine properties of the tagged neutron beam by identical 
cuts apphed to event streams (l)-(3), so that associated 
systematic uncertainties would cancel in the yield ratios 
from which the absolute np scattering cross section is ex- 
tracted. Among these common cuts, described in more 
detail below, were ones to remove BD noise contribu- 
tions correlated among the four quadrants of the tagger, 
to identify the recoil particles detected in the tagger, and 
to divide the tagged neutron events into subsamples for 
subsequent analysis. Additional cuts defined a fiducial 
range for the tagged neutron's predicted transverse co- 
ordinates at the secondary target (| xtag |<9.5 cm and 
I ytag |<9.5 cm), and selected prompt tagger two-particle 
coincidences {\tpi — tp2 — 30 ns| < 70 ns, where tpi (tp2) 
is the arrival time of the recoil proton with the larger 
(smaller) DSSD energy deposition). Software cuts ap- 
plied to event stream (2) alone to identify free np scat- 
tering events were kept to a minimum in order to avoid 
complicated systematic errors. We relied instead on the 
accuracy of the background subtractions, which could be 
verified to high precision. Before application of cuts, ad- 
ditional MWPC requirements were added in software to 
amplify the hardware definitions of the various triggers. 
Thus, at least one hit in the x-plane and at least one hit in 
the j/-plane were required for each of the three MWPC's 
for events from stream 2. 



1. Particle identification 

The correlation of DSSD vs. BD energy depositions 
was used to select two basic event classes for analysis 
of each of the three tagged neutron event streams: (a) 
"2-stop" events, where both protons associated with the 
neutron stopped inside the DSSD (either the same or 
different quadrants of the tagger); and (b) "1-punch" 
events, where one of the protons stopped inside a DSSD 
and the other punched through to the BD in a differ- 
ent quadrant and stopped there. These two classes, as 
discussed further below, differ significantly in neutron en- 
ergy (En) and position (xtag) profiles, allowing an impor- 
tant crosscheck on the accuracy of the tagging technique 
by comparing np cross sections extracted independently 
from each class. The 2-stop events were further sub- 
divided according to whether the higher of the two re- 
coil proton DSSD energy depositions (Epi) was below or 
above 5.0 MeV. Protons with Epi > 5.0 MeV range out 
near the exit of the DSSD, and hence possibly in dead 
layers at the back of the DSSD or front of the BD, making 
this event class subject to somewhat greater ambiguity 
regarding complete recoil proton energy reconstruction 
and accuracy of the predicted Xtag value for the tagged 



neutron. Events where both protons punched through 
to the BD's, or where either punched through the BD 
itself, were not included in the analysis because they cor- 
responded to En below the range of interest. 
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FIG. 3: Particle identification plot for one DSSD-BD combi- 
nation for (a) tagged protons and (b) tagged neutrons. 

Figure 131 shows raw spectra for both tagged neutrons 
(event stream 1 in frame (b)) and tagged protons (event 
stream 4 in frame (a)) of the energy deposited in a typi- 
cal DSSD quadrant vs. that in the companion BD when 
the latter is non-zero. The tagged proton events exhibit 
clear recoil proton (lower) and deuteron (upper) particle 
identification loci, while only the proton locus remains 
for tagged neutron events. The loci bend backwards 
when the detected particle begins to punch through the 
BD. The two-dimensional gate (dark boundary) shown 
in each frame was used to select recoil protons that en- 
ter and stop inside the BD, e.g.^ to identify the 1-punch 
tagged neutron events. Note that the most intense re- 
gion along the proton locus, corresponding to deuteron 
breakup events with an energetic large-angle proton, is 
thereby eliminated. So are events lying off the proton 
locus, where the backing detector response may be cor- 
rupted by noise or pileup. 

Figure 21 shows the reconstructed En and Xtag distri- 
butions for the tagged neutrons in the 2-stop (for all val- 
ues of Epi) vs. 1-punch samples. While the two sam- 
ples yield overlapping distributions, it is clear that the 

1- punch events correspond on average to lower-energy 
neutrons at larger production angles (preferentially pop- 
ulating the beam-right side of the secondary target). 

2. Correlated noise in the BD 

Special care was taken in the definitions of 1-punch and 

2- stop events to minimize effects of substantial detector 
noise picked up by the large-capacitance BD's. An im- 
portant source of this noise was discovered to arise from 
the initiation of pulse height information readout on the 
adjacent DSSD front-end electronics 0|. The induced 
noise was strongly correlated among the four BD's, as 
illustrated in Fig.|21 This figure reveals two uncorrelated 
bands parallel to the x and y axes, due to 1-punch events 
in one of the quadrants and low pulse-height noise in the 
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FIG. 4: The reconstructed energy (a) and horizontal transverse coordinate (b) of tagged neutrons at the secondary target for 
1-punch events, 2-stop events and their sum (sohd black line). 



other (the pedestals for each BD appear in ADC channel 
«10). But one also observes a strong diagonal band in- 
dicative of noise correlations between the two quadrants. 
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BD Quadrant 4 ADC Channel # 

FIG. 5: The ADC outputs of two BD's plotted against one 
another, revealing regions of 1-punch tagged neutron events, 
but also a region dominated by correlated noise in both quad- 
rants. Events along the x- and y-axes, where one BD sig- 
nal falls in the pedestal region, are also mostly valid 1-punch 
events. 

Since the noise correlation pattern extends beyond a 
reasonable software threshold, it was necessary to use 
a two-dimensional gate, such as that shown in Fig. 
to bound the noise correlation region. Candidates for 
valid 1-punch events were then required to: (1) surpass a 
threshold ADC channel (sal5) on at least one BD; (2) fall 
outside the noise correlation gates for all BD pairs; (3) 
not surpass the BD noise peak (ADC «70) in more than 



one BD; and (4) fall within the PID gate in Fig.|31for the 

appropriate quadrant. These conditions and the comple- 
mentary ones required for 2-stop events reduced the flux 
of tagged neutrons considered for subsequent analysis, 
by removing events with BD pulse height ambiguities, 
but since they were applied equally to all tagged neutron 
event streams, they did not introduce systematic errors 
in the np cross section extraction. 

B. Corrupted events subtraction 

An event misidentification mechanism discovered dur- 
ing the data analysis was attributed to an electronics 
malfunction in the gating or clearing circuit for the elec- 
tronics module that was used to digitize the pulse height 
information for all four BD's. The effect of the mal- 
function was to zero out valid BD energy signals for 
a randomly selected fraction of punch-through events. 
The effect was seen clearly, for example, in the pd elas- 
tic scattering events in stream 4, where a software gate 
placed on the two-body kinematic correlation between re- 
coil deuteron DSSD energy deposition and forward pro- 
ton angle could be used to select events in which the 
deuteron must have stopped in the BD. Roughly 3/4 of 
these events showed the anticipated BD pulse height, but 
1/4 had Ebd = 0. In the case of tagged neutrons, the 
corrupted events were misidentified as 2-stop events and 
gave systematically incorrect predictions of the tagged 
neutron trajectory, since some recoil proton energy was 
lost. However, the availability of full information for the 
surviving punch-through samples allowed us to emulate 
the effect and subtract the corrupted events accurately. 

The corrupted events were easily distinguished in event 
stream 2 by using the MWPC tracking information. Fig- 
ure El shows the correlation for event stream 2 between 
Epi in the tagger and Xtrack-xtag, where Xtrack denotes 
the transverse coordinate of the detected proton from 
np scattering at the secondary target, as reconstructed 
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from the MWPC hits. The majority of events have 
Xtrack — Xtag ~ 0, independent of Epi, as expected when 
both the tagging and tracking are accurate. The cor- 
rupted events populate the "tail" to the left of the most 
intense band, with Xtag exceeding Xtrack by an amount 
that is strongly correlated with the recoil proton energy 
deposition in the DSSD: the lower Epi, the larger is the 
lost Ebd and the consequent error in xtag- While the cor- 
rupted events could be eliminated from event stream 2 by 
a software gate within Fig. they could not have been 
similarly eliminated from event streams 1 and 3, where 
there was no forward proton to track. Hence, it was es- 
sential to find a way to subtract these corrupted events 
reliably and consistently from all three tagged neutron 
event streams. 
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FIG. 6: The correlation for apparent 2-stop np scattering can- 
didates between the higher of the DSSD energy depositions 
for the two recoil protons in the tagger and the difference in 
predicted x coordinates at the secondary target from neutron 
tagging vs. forward proton ray-tracing. The long correlated 
tail contains corrupted punch-through events for which elec- 
tronic loss of backing detector energy information has led to 
misidentification of the event and large systematic errors in 
the tagging. 



The corrupted events were simulated using all recorded 
punch-through events that survived with their BD energy 
information intact, by reanalyzing these events after ar- 
tificially setting Ebd = in the software before tagging 
reconstruction. The distribution shapes of the tail events 
in Fig. 13 with respect to all variables were accurately re- 
produced when this simulation was based on all events 
in Fig.|2Ib), both inside and outside the two-dimensional 
gate drawn, and also including events where both recoil 
protons punched through their DSSD's. In order to de- 
termine the fraction of these punch-through events that 
was affected by the electronics malfunction, we relied on 



a comparison of the subsamples of our simulated events 
and of the apparent 2-stop events that had valid BD tim- 
ing information despite having Ebd = 0. Because the 
BD noise problems necessitated high thresholds to gen- 
erate timing signals, these subsamples populate mostly 
the far tail in Fig. corresponding to Epi < 3.5 MeV 
(thus, to relatively large BD analog signals). The cor- 
rupted fraction of punch-through events was in this way 
determined independently for each of the three tagged 
neutron event streams, and found to be identical for the 
three, within the statistical precision (typically « 1%) 
available in matching simulated and recorded corrupted 
event subsamples. The fraction varied slightly with time 
during the production run, but averaged 23%. The suc- 
cess of this simulation and the persistence of the cor- 
rupted fraction across (both tagged neutron and tagged 
proton) event streams provide strong support for our as- 
sumption that the malfunction affected a random sample 
of punch-through events. Because the fractional loss of 1- 
punch events to this corruption was independent of event 
stream, there was no residual systematic effect on the ex- 
tracted 1-punch cross sections, but rather only a slight 
loss of statistical precision. 
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FIG. 7: The subtraction of the simulated corrupted events 
removes the tail from the 2-stop event sample. The simulated 
spectrum has been normalized by matching to the 2-stop sub- 
sample characterized by backing detector information with 
valid times but zero energy. 



Figured shows a comparison of the x track — Xtag spec- 
trum for all 2-stop events in stream 2 with the simu- 
lated corrupted sample, normalized as described above 
via the subsample with valid BD timing signals. The sub- 
traction eliminates essentially completely the corrupted 



events with Epi < 5.0 MeV, or xtra 



Xtac 



< — 5 cm, 



leaving a reasonably symmetric small background (dis- 
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cussed further in Sec. IV AS)) at \xtrack — Xtag\ > 5 cm. 
We therefore assume that the subtraction is similarly suc- 
cessful for event streams 1 and 3, where we have no track- 
ing information to compare, and associate a systematic 
error for the subtraction (see Sec. IVB 2|l that reflects 
only the uncertainty in the normalization scheme for the 
simulated corrupted events. 

For Epi > 5.0 MeV, there is a remaining tail of small 
extent in the subtracted Xtrack ~ Xtag spectrum in Fig.[7| 
that arises not from the electronics malfunction, but 
rather from recoil protons that barely punch through the 
DSSD, while depositing insufficient energy in the BD to 
be distinguished from noise. Because of these events, we 
have separately analyzed the 2-stop samples with Epi < 
5.0 MeV and Epi >5.0 MeV. For the latter sample, after 
subtracting simulated corrupted events, we used a two- 
dimensional software gate on Fig.^jto eliminate events in 
stream 2 that had potentially distorted xtag information, 
thereby rejecting 18% of 2-stop Epi >5.0 MeV events 
(as opposed to the 23% of all 2-stop events in stream 2 
that were affected by the corruption). The yields of 2- 
stop Epi >5.0 MeV events in streams 1 and 3 were then 
scaled down by the same 18%, to remove the remaining 
events of questionable 2-stop pedigree. 

The small peak at Xtrack — Xtag ~ in the simulated 
background in Fig.[7|indicates that a small fraction of the 
punch-through event sample used in the simulation really 
corresponds to true 2-stop events that were misidcntified 
by virtue of BD noise that evaded the noise cuts discussed 
in the preceding subsection. Subtracting this small frac- 
tion of valid 2-stop events along with the simulated cor- 
rupted events has the effect of reducing the 2-stop tagged 
neutron yield by « 3% in all three event streams, with no 
significant consequence for the absolute np cross sections 
extracted from the 2-stop sample. 



C. Background subtraction 

The background events for this experiment came 
mostly from np quasifree scattering off carbon nuclei in 
the CH2 target. However, there were also some promi- 
nent sources displaced from the secondary target, includ- 
ing: (1) protons coming directly from the gas jet produc- 
tion target, or from the exit flange on the Cooler beam 6° 
magnet chamber (see Fig.|5}, that evaded the veto scintil- 
lators due either to their imperfect coverage or electronic 
inefficiencies; (2) np scattering events induced either on 
the downstream scintillator face or the Lucite light guide 
for the SUV, yielding pulse heights below that veto detec- 
tor's threshold; and (3) quasifree np scattering induced 
on the vertically narrow (but longitudinally thick) alu- 
minum frame used to support the secondary target. By 
frequently interchanging the CH2 target with a graphite 
target closely matched in transverse dimensions and in 
areal density of carbon nuclei, we were able to subtract 
the backgrounds from all sources simultaneously. The 
relative normalization of the CH2 and C runs was deter- 



mined from the pd elastic scattering yield from the GJT, 
as recorded in event stream 4 17J. The background sub- 
traction was determined to be sufflciently reliable that 
we could avoid imposing many kinematic cuts, with po- 
tentially significant systematic ambiguities, to define free 
np scattering events. 

The accuracy of the background subtraction can be 
judged, for example, from Fig. |S1 which presents CH2 
and C spectra, and their difference, with respect to ytag 
(the vertical impact position of the neutron on the sec- 
ondary target, as reconstructed from the tagger) and AE 
scintillator pulse height within a narrow np scattering 
angle range. These two particular variables have been 
chosen for display in the figure because the CH2 spec- 
tra show prominent background features associated both 
with quasifree scattering (the long high pulse height tail 
in AE) and with other sources (the ytag = — 11 cm peak 
from the aluminum support frame). Both sources are 
precisely eliminated by the background subtraction. In- 
deed, upper limits on the surviving remnants of these 
features allow us (as described in Sec. IVB l|l to reduce 
the systematic uncertainty for the background subtrac- 
tion even below the level given by the precision of the 
measured C/CH2 target thickness ratio (±0.6%). 



D. Free scattering cuts 

Software conditions imposed only on event stream 2 to 
distinguish free-scattering from background events have 
the potential to remove free-scattering yield in sometimes 
subtle ways. They were thus used in the analysis only 
when they could substantially reduce the statistical un- 
certainties (i.e., by suppressing background to be sub- 
tracted) without introducing significant systematic un- 
certainties in correcting for the free-scattering losses, or 
when such losses were judged to be inevitable to remove 
ambiguities in the analysis. The accuracy of the C back- 
ground subtraction provided a reliable method to judge 
the extent of any free-scattering event removal. 

The most effective such cut applied was placed on the 
correlation of forward proton energy loss in the AE scin- 
tillator with the laboratory angle of the proton trajec- 
tory. The applied two-dimensional software gate is su- 
perimposed on the observed distribution of events fol- 
lowing CH2 - C subtraction in Fig. O This distribution 
reveals that very few free-scattering events were removed 
by this gate, but it is clear from the long tail seen in 
the projected unsubtracted spectrum for one angle bin 
in Fig. Oc) that a substantial number of quasifree back- 
ground events, leading to lower-energy outgoing protons, 
were successfully removed. 

In contrast, we did not apply a comparable cut on 
the energy deposition of the forward proton in the rear 
hodoscope, where it generally stopped, despite an appre- 
ciable difference in the distributions of hodoscope energy 
between free and quasifree events. The reason for avoid- 
ing this cut is illustrated in Fig. ^1 the free-scattering 
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FIG. 8: Distribution of np scattering event candidates with respect to ytag and to AE scintillator pulse height for the CH2 and 
C targets separately, and for the difference between them. 
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FIG. 9: The distribution of np scattering candidate events 
after subtraction of C from CH2 data with respect to AE 
ADC and proton lab angle. The red lines show the boundaries 
of the gate applied to event stream 2 to select free scattering 
events. 
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FIG. 10: The distribution of np free scattering events in the 
scattering angle bin 9p'^ = 3 — 6° with respect to forward 
proton energy deposition in the rear hodoscope. The curve 
represents a fit with a Gaussian plus exponential tail. The 
tail represents valid free scattering events where the proton 
undergoes a nuclear reaction in the stopping scintillator. 



spectrum revealed by the C subtraction exhibits a quite 
substantial low-energy reaction tail in addition to the 
well-defined full-energy peak. An unacceptably large sys- 
tematic error would have been introduced by the need to 
correct for loss of these reaction tail events, if we had im- 
posed a cut on hodoscope energy to suppress background. 

A benign cut imposed for slight background reduction 
placed an upper threshold on the value obtained for a 
linear track fit to the reconstructed MWPC space points. 
The CH2 - C subtraction indicates that only (0.2 ±0.1)% 
of free-scattering events were removed by this condition. 
More serious (6.3% of total CH2 - C yield), but unavoid- 
able, losses were introduced by cuts confining the tag- 
ging and tracking information in event stream 2 to agree 



within \xtrack - Xtagl < 2.5 Cm and \ytrack - ytagl < 2.0 

cm. These limits correspond to ±3cr of the narrow Gaus- 
sian resolution function that dominates these distribu- 
tions in the CH2 - C spectra. Nonetheless, the cut elim- 
inates events in long distribution tails that are affected 
either by tagging errors or sequential reactions of the 
tagged neutron, which introduce serious ambiguities in 
interpretation. This cut, and its consequence for system- 
atic uncertainties, will be discussed further in Sec. IV A'Sl 
Finally, it is worth mentioning one additional cut that 
we chose not to impose. The transverse np vertex coor- 
dinates are, in fact, determined by the tagging and track- 
ing with considerably better resolution than implied by 
the CT w 7 — 8 mm value mentioned in the preceding 
paragraph 0. This latter value is dominated by the 
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thickness of the secondary target, simply reflecting the 
uncertainty in precise longitudinal origin of the np scat- 
tering vertex. Much better information is, in principle, 
available by locating the vertex in three dimensions at 
the point of closest approach of the neutron trajectory 
reconstructed from the tagger and the proton trajectory 
reconstructed from the MWPC's. Distributions of such 
reconstructed secondary vertex coordinates 0| permit 
tagged neutron radiography of the background sources 
displaced from the CH2 target. However, any cuts to re- 
move such background sources in this manner would be 
affected by the strong dependence of the reconstructed 
vertex resolution on the proton scattering angle (ver- 
tex information clearly deteriorates as the neutron and 
proton trajectories approach coUinearity) . We decided 
to rely completely on the C subtraction to remove such 
other sources of background, in order to avoid consequent 
angle-dependent free-scattering event losses. 



E. Acceptance 

The lab-frame proton scattering angle is deter- 
mined for each analyzed event as the opening angle be- 
tween the neutron trajectory reconstructed from the tag- 
ger and the forward proton trajectory reconstructed from 
the MWPC's. The geometric acceptance of the forward 
detector array for np scattering events is a function of 
both and the coordinates of the scattering vertex at 
the secondary target. Because the distribution of scat- 
tering vertex coordinates, especially of Xtag (see Fig. 0J, 
differed among the three analyzed data subsamples (1- 
punch, 2-stop with Epi < 5.0 MeV, and 2-stop with 
Epi > 5.0 MeV), the acceptance had to be evaluated 
separately for each subsample. This was done by com- 
paring simulated to measured distributions of events with 
respect to azimuthal angle (j)^'^ within each 6^'^ bin, for 
each data subsample. 

The simulations were constrained to reproduce the 
measured distributions of the longitudinal production 
vertex coordinate of the neutron within the GJT (com- 
mon to all three data subsamples) and its transverse co- 
ordinates on the secondary target (separately for each 
subsample). Within these distributions, coordinates were 
generated randomly for each event, as were also Op'^ (in 
the range 0-75°) and 0*'^ (over the full azimuthal range). 
Generated outgoing proton trajectories were then ac- 
cepted if they would yield signals above the hodoscope 
pulse height threshold (required in trigger) and in all 
three MWPC's (required in the data analysis). Forward 
detector location parameters were tuned slightly from 
their measured values to optimize the fit of the simu- 
lated to the measured ^p'^ distributions for all 9p'^ bins 
and for 1-punch and 2-stop samples simultaneously. 

The high quality of the fits obtained is illustrated in 
Fig. ^2 for the 1-punch (a) and 2-stop (b, summed over 
all Epi) samples for a single large angle bin, 6*^^=42-45°, 
where the observed azimuthal distributions display con- 
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FIG. 11: Comparison of the measured (solid line with sizable 
statistical fluctuations) and simulated (dotted line) distribu- 
tions of free (CH2-C) np scattering events in the angle bin 
9p'^ — 42 — 45° with respect to proton azimuthal scattering 
angle (fip'^, for the 1-punch (a) and 2-stop (b) data samples. 
Forward detector geometry parameters, plus a single overall 
normalization parameter per angle bin and data sample, have 
been adjusted to optimize the fit simultaneously for all angle 
bins and both data samples. 



siderable structure. The structure reflects the rectangu- 
lar shape of the hodoscope and large MWPC, projected 
onto 9 — (f> space: e.g.^ the four peaks observed correspond 
to the four detector corners. The small changes in dis- 
tribution between the 1-punch and 2-stop samples - e.g., 
in the relative heights of the peaks and in the extent of 
the dips near ipp'^ — 0° (beam left side) and 180° (beam 
right) - arise from the shift in xtag profiles seen in Fig. 21 
These features are all reproduced very well by the sim- 
ulations. For 9p'^ < 24°, the measured and simulated (j) 
distributions are essentially uniform over 27r, indicating 
full acceptance. Figure El shows the simulated accep- 
tance for the 1-punch data sample as a function of 9p'^. 
The 0.2% shortfall from full acceptance near 0° reflects 
protons incident normally on the small cracks between 
adjacent hodoscope elements. Results presented in the 
next section are limited to the angle range for which the 
acceptance is at least 50%; at larger proton angles the 
uncertainty in acceptance grows rapidly. 



IV. RESULTS 

The absolute differential cross section for np back- 
scattering was extracted independently for three data 
samples - 1-punch, 2-stop with Epi <5 MeV, and 2-stop 
with Epi >5 MeV - from the yields of event streams 1, 
2 and 3 as follows: 

.da ^ N,{9;^)Uc^ 

^dn'""" iNi + N2 + N,)tH\dcos{9^p-)\a49;-)' ^> 

where A^2(^p'^) represents the number of free-scattering 
events from stream 2 within a given reconstructed pro- 
ton angle bin, surviving all relevant cuts and background 
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TABLE I: Correction factors and systematic uncertainties in 
correction factors for the np cross sections. 



FIG. 12: The simulated acceptance of the 1-punch data sam- 
ple. The inset shows a greatly magnified vertical scale for the 
most forward proton scattering angles. 



subtractions; A^i, iV2 and N3 in the denominator rep- 
resent analogous tagged-neutron yields from the mutu- 
ally exclusive event streams 1 (corrected for prescaling), 
2 (angle-integrated) and 3; the Ci represent small cor- 
rections, summarized in Table ^ with details in section 
IVI for various inefficiencies, tagged neutron losses or 
backgrounds, and software cut and dead time differences 
among event streams; tn = (1.988 ± 0.008) x 10^3 H 
atoms/cm^ for the CH2 target; and is the azimuthal 
acceptance determined from simulations for the given an- 
gle bin. The data were analyzed in 1 MeV wide slices of 
reconstructed neutron energy from 185 to 197 MeV and 
an effective cross section extracted at the mean neutron 
energy of 194.0±0.15 MeV. For this purpose, a small (al- 
ways < 1%) cross section correction was made for the 
deviation of each analyzed slice from the mean energy, 
using the theoretical energy- and angle-dependence cal- 
culated with the Nijmegen PWA93 solution j^. 

The cross sections for the three data subsamples, with 
their independently determined absolute scales, are mu- 
tually consistent in both magnitude and angular shape, 
within statistical uncertainties, as revealed by the com- 
parisons in Fig. ^1 The figure shows the relative differ- 
ence, {{"^j^^ sampleA {'^]^\sampleB) / {"^^q) sampleB : between 

pairs of cross sections for the three data samples. The 
reduced value for the comparison of each pair of sam- 
ples is indicated in the legend of Fig. El This compari- 
son supports the reliability of the experiment and analy- 
sis, because these samples come from complementary re- 
gions of the tagged beam spatial and energy profiles (see 
Fig. 0J and are subject to somewhat different system- 
atic error concerns. We view the agreement in absolute 
cross section scale as particularly significant demonstra- 



Source 


Correction Factor (ci~ 


Uncertainty in Ci 


Accid. tagger coinc. 


1.0003 


< ± 0.001 


Non-D2 tagger 


1.0067 (2-stop); 


± 0.002 


background 


1.0044 (1-punch) 




n pos'n unc. on CH2 


1.0000 


± 0.001 


n atten'n before CH2 


1.005 


± 0.0025 


Sequential react'ns 


1.063 


± 0.010 


& a^ta<,(n) errors 






C bkgd. subtraction 


1.0000 


± 0.004 


H in C target 


1.000 


± 0.004 


Corrupted event 


1.000 


< ±0.001 


subtraction 






Software cut losses 


1.010 


± 0.005 


Reaction tail losses 


1.004 


± 0.002 


Neutron polarization 


Angle-dependent: 


± 0.001 


effects 


< 1.0012 (1-punch) 










CH2 tgt. thickness 


1.0000 


± 0.004 


np scattering 
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acceptance 




±0.017 (90°) 


MWPC inefficiency 


1.017 
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Trigger inefficiency 
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cos2(6l^^^) 


X cos2(e^^^)] 


Dead time diffs. 


0.991 


± 0.005 


Scattering angle 


1.000 


angle-dependent , 


errors 




< ±0.004 


Net, typical 


fa 1.10 


^ ±0.016 



tions of the accuracy of the neutron profiles reconstructed 
from tagging and of the subtraction procedure applied 
to remove corrupted events from the 2-stop sample (see 
Sec. III.B). Cross sections extracted for different time pe- 
riods within the production runs, and with different sets 
of cuts, are also consistent within uncertainties. 
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FIG. 13: The fractional differences between the absolute dif- 
ferential cross sections extracted for different analyzed data 
subsamples. The plotted error bars take into account only 
the independent statistical (including those from background 
subtractions) uncertainties for the three samples. Slightly dif- 
ferent correction factors Ci were applied to the cross sections 
for different samples, as indicated in Table Q before the com- 
parison was made. 
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FIG. 14: Absolute differential cross section from the present 
experiment, compared with data from Ref. i3] and with PWA 
calculations at two relevant energies. The error bars on the 
present results are statistical (including background subtrac- 
tion), while the shaded band represents all systematic uncer- 
tainties, including those in the overall normalization. 

The results, averaged over all three data samples, are 
compared in Fig. ll4l with previous experimental results at 
162 MeV (2| and with the Nijmegen partial wave analysis 
(PWA93) at the two relevant energies Q . The measured 
points are plotted at the yield-weighted centroid angle 
of each analyzed bin. The comparison of the present 
results with previous experiments and with partial wave 
analyses will be discussed in detail in Sec. VI, after first 
describing the nature and evaluation procedure for each 
of the systematic uncertainties included in Table ^ 

V. SYSTEMATIC ERRORS 

Most of the individual correction factors ci applied to 
the extracted cross sections, and their associated system- 
atic uncertainties listed in Tabled have been evaluated 
via complementary analyses of the data. In this section 
we briefly describe the procedures used and error esti- 
mates for each, being careful to distinguish uncertainties 
that affect only the overall cross section normalization 
from those with appreciable angle-dependence. In the 
latter cases, we also characterize the degree of correla- 
tion among the uncertainties at different angles, to fa- 
cilitate inclusion of the uncorrelated systematic errors in 
PWA's including the present data. For purposes of logi- 
cal flow, we organize the discussion into four categories: 
(1) tagged neutron flux uncertainties; (2) np backscatter- 



ing yield uncertainties; (3) target thickness, acceptance 
and efficiency uncertainties; and (4) errors in determining 
the kinematic variables. 



A. Tagged neutron flux uncertainties 

The sources below contribute to uncertainties in 
extracting the angle-integrated yields Ni^2,3 in Eq- Q, 
dominated by the non-interacting tagged neutrons in 
event stream 1. All of the issues discussed in this 
subsection give rise to overall (angle-independent) 
normalization errors in the differential cross sections. 



1. Accidental tagger coincidences 

Accidental coincidences between two uncorrelated par- 
ticles detected in the tagger contribute slightly to the 
apparent tagged neutron flux on the secondary target, 
leading to an underestimate of the cross section. The 
accepted events in all three event streams passed a cut 
on the time difference At = (tpi — tp2) between the two 
tagger hits, as indicated in Fig. E| The correction fac- 
tor was determined from the ratio of events in stream 
1 that passed all other cuts defining the tagged neutron 
beam but fell within one of two displaced time windows, 
I At + 110 ns |< 70 ns and | At - 170 ns |< 70 ns, to the 
yield in the prompt coincidence window | At— 30 ns |< 70 
ns. The resulting correction factor is ci = 1.0003, with 
an uncertainty < ±0.001, showing that accidental coin- 
cidences were a minor issue for the experiment. 



2. Tagger background from non-D2 sources 

Additional possible background contributions to the 
tagged neutron flux could arise from real (correlated) 
two-particle coincidences in the tagger, generated by pro- 
ton beam interactions with nuclei heavier than deuterium 
in material displaced from the GJT. This possibility was 
checked via runs where H2 was substituted for the D2 
in the GJT, to induce similar beam "heating" without 
any real possibility of tagged neutron production (since 
the proton beam energy was far below pion produc- 
tion threshold for the p+p system). A correction fac- 
tor C2 = 1.0044 ± 0.002 (1.0067 ± 0.002) for 1-punch (2- 
stop) events was determined from the ratio of accidental- 
subtracted tags satisfying the tagged neutron conditions 
with the H2 vs. D2 production targets. The statistical 
uncertainties in these ratios were considerably smaller 
than ±0.002; the quoted uncertainty is intended to allow 
for the possibility of slight systematic differences in beam 
heating, hence in the rate of interactions with displaced 
material, between the two GJT gases. 
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FIG. 15: The distribution of DSSD arrival time difference 
between the two recoil protons for tagged neutron events, 
shown for all 2-stop events (following subtraction of the cor- 
rupted events as per Sec. If If Bt and for 2-stop events where 
the lower of the two DSSD energy depositions (Ep^) exceeds 
2.0 MeV. The vertical lines indicate the prompt coincidence 
gate (central region between the inner two lines) imposed in 
the analysis, plus two displaced gates used to assess acciden- 
tal coincidence background in event stream 1. The long tail 
seen for events with Ep2 < 2.0 MeV, arising from detector 
noise and imperfect software corrections for time walk near 
the front-end discriminator threshold, leads to an overesti- 
mate of the accidental coincidence yield, but the correction 
and uncertainty still remain quite small. 



3. Impact position uncertainty on the CH2 target 

This is the first of several error sources we consider that 
arise when a properly tagged neutron does not reach the 
CH2 (or C) secondary target, or reaches it at a signifi- 
cantly different position than expected from the tagging. 
Because of the finite (several mm) impact position reso- 
lution from the tagger, some tagged neutrons predicted 
to hit the secondary target may actually miss it, while 
some predicted to miss the target may hit it. Especially 
near the target edges, where the yield of np scattering 
events drops rapidly and nonlinearly as a function of im- 
pact position, this resolution smearing can affect the ex- 
tracted cross sections. In practice, however, we observe 
no statistically significant difference in cross section nor- 
malization between the independent event samples from 
the bulk of the target (| xtag |< 9.0 cm and | ytag |< 9.0 
cm) and from a 5.0-mm wide strip (9.0 <| Xtag |< 9.5 
cm or 9.0 <| ytag |< 9.5 cm) surrounding this core. 
From this comparison and the fraction of all events aris- 
ing near the target edges, we infer a correction factor 



C3 = 1.000 ±0.001. 



4- Neutron attenuation before the CH2 target 

Some tagged neutrons fail to hit the secondary target, 
leading to an underestimate of the extracted np cross sec- 
tion, as a result of interactions they undergo upstream 
of that target. Approximately 3.5% of 200 MeV neu- 
trons will undergo an inelastic reaction of some sort in 
the upstream material ,25i |. which is dominated by the 
0.29 cm thick stainless steel vacuum window at the exit 
of the Cooler's 6° magnet vacuum chamber, the 0.64 cm 
thick LUV plastic scintillator, and the 0.64 cm thick SUV 
plastic scintillator (the first two of these traversed at an 
incidence angle ~ 14°). However, many of these "pre- 
scattering" neutrons give rise to charged products that 
get vetoed by LUV or SUV (and hence do not contribute 
to the tagged flux) or are removed by the C subtrac- 
tion as apparent np scattering events from an upstream 
source. Others yield an energetic neutron or proton, not 
strongly deflected from the original tagged neutron tra- 
jectory, that still strikes the nearby secondary target with 
a chance to induce a tertiary interaction there, and so 
might still be considered as part of the incident flux. 

We may judge the rate of such tertiary interactions 
from events where a forward proton emerges from the 
secondary target at a transverse location {xtrack.y track) 
substantially different (by much more than the tagging 
position resolution effect considered in Sec. IV A3|) from 
the predicted impact position of the tagged neutron 
{xtag,ytag)- Such tertiary interactions introduce their 
own problems in the analysis, to be addressed separately 
in the next subsection. Here, we study them in order to 
place limits on the probability of larger-angle upstream 
scattering, which yields no chance of a tertiary scattering. 
In Fig.^1 we show the difference spectra for Xtrack — xtag 
and y track — ytag for 1-punch events (to eliminate ambigu- 
ities from the corrupted 2-stop events seen in Fig.O after 
CH2 - C subtraction (to eliminate ambiguities from np 
scattering induced on material displaced from the sec- 
ondary target). The narrow Gaussian resolution peaks 
sit atop a broad background that has important contri- 
butions from these tertiary interactions. 

By fitting the background in Fig. 1161 with a broad 
Gaussian, and assuming that the probability of initiat- 
ing a scattering in the secondary target is roughly the 
same for the pre-scattered neutrons as for the bulk of the 
tagged neutrons, we estimate that 0.5% of tagged neu- 
trons may be pre-scattered through a sufficiently large 
angle to cause a transverse displacement greater than 10 
cm [i.e., half the target width or height) on the secondary 
target. We use this estimate to infer a correction factor 
C4 = 1.005 ± 0.0025 for tagged neutron pre-scattering 
flux losses before the CH2 target. The ±50% uncer- 
tainty we assign to (c4 — 1) is intended to account for 
non-prescattering origins of the background in Fig. 1161 
(e.g., tagging errors or sequential reactions within the 
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CH2 or tertiary scattering of a forward proton in mate- 
rial following the CH2 target) and for possible upstream 
neutron interactions that elude the above analysis. 




FIG. 16: The diflerence distributions in x (left) and y (right) 
coordinates determined from tracking of the forward proton 
vs. reconstruction of the tagged neutron, for the 1-punch data 
sample after CH2 - C subtraction. The vertical lines indicate 
the location of software gates used to remove events that may 
be complicated by sequential reactions or tagging errors. 



5. Sequential reactions in the secondary target or upstream 
material 

Here we deal explicitly with the events contributing 
to the broad backgrounds in Fig. ^| and in the analo- 
gous distributions for 2-stop events. From the behav- 
ior of these events and the differential cross sections we 
extract specifically from them, we judge them to corre- 
spond primarily to valid free np scattering in CH2 either 
following an earlier interaction of the tagged neutron or 
preceding a later interaction of the forward proton. The 
two interactions will in some cases both have taken place 
within the CH2 target. Some of the background may 
also arise from tagging errors associated with less than 
complete energy collection for the recoil protons in the 
tagger. Regardless of their detailed source, such events 
are distorted because we mismeasure the scattering angle 
and possibly the incident neutron energy for the free np 
scattering. 

The least biased way to handle these events is to elim- 
inate them from the analyzed sample. This is simple 
enough to do for each angle bin in event stream 2, via 
the software cuts requiring \x track — Xtag\ < Sax and 
\ytrack - ytagl < ScTy, whcrc (Ta; ~ 0.8 Cm and cry « 0.7 cm 
are the widths of the narrow Gaussian peaks in Fig. 1161 
However, we have no access to analogous cuts for event 
streams 1 and 3, where there is no MWPC information. 
Hence, we must correct the extracted cross section for our 
inclusion of tagged neutron flux that is associated with 
such eliminated sequential reaction events. The correc- 
tion factor C5 = 1.063 ± 0.010 is the largest single cor- 
rection and systematic uncertainty we apply. The value 
assumes that the 6.3% of event stream 2 (CH2 - C) yield 
(averaged over 1-punch and 2-stop samples and over scat- 
tering angle) removed by these sequential reaction cuts 



arises from 6.3% of the tagged neutron flux. (The ac- 
tual correction factors applied differ for the three data 
samples, reflecting differences in the fraction of events 
removed by these cuts.) The uncertainty allows for er- 
rors in this assumption, for example, because it neglects 
the energy-dependence of the np scattering probability 
in the CH2 target. 



B. Uncertainties in absolute np backscattering 
yields 

Analysis issues in the extraction of the free-scattering 
yield N2iO) needed in Eq. □ can lead, in principle, to 
angle-dependent errors. We thus specify for each case 
below whether the estimated uncertainty should be con- 
sidered as angle-dependent and as uncorrelated from an- 
gle bin to bin. 



1. Uncertainties in background subtraction via the C target 

As described in Sec. 1111 CI we relied heavily on the 
CH2 - C subtraction to remove simultaneously back- 
grounds due to quasifree scattering from carbon nuclei 
in the secondary target and to reactions induced on dis- 
placed sources. The precision of the subtraction depends 
on that of our knowledge of the relative target thick- 
nesses and integrated neutron flux exposures for the CH2 
vs. C runs, and on the stability of beam conditions be- 
tween the two sets of runs. The relative normalization, 
taken from cleanly (kinematically) identified pd scatter- 
ing yields measured simultaneously, is determined with 
quite high statistical precision but could, in principle, 
deviate systematically from the more relevant ratio of 
tagged neutron yields. The overall precision of the rela- 
tive normalization was judged from the extent to which 
scattering events from the aluminum target platform (see 
Fig. IHIa)) were successfully removed by the C subtrac- 
tion. We concentrate first on this background source 
because its yield is not sensitive to the CH2/C target 
thickness ratio. 

The reconstructed ytag distributions in the vicinity of 
the aluminum platform peak, for both CH2 and C tar- 
gets (see Fig. |Sl , could be well reproduced by the sum 
of a Fermi distribution and a polynomial, to represent 
the bottom target edge, and a Gaussian to represent 
the aluminum peak. The fit for the CH2 is shown in 
Fig. lTTT a'). An analogous fit was then carried out for the 
C-subtracted spectrum in Fig. |Htb) , fixing the positions 
and widths of the Gaussian and Fermi-function contribu- 
tions to their common values for CH2 and C. The ratio 
of events in the Gaussian peak after subtraction to that 
before subtraction is (1.9 ± 0.54) x 10"'^, providing one 
measure of the accuracy of the background subtraction. 

An independent measure was provided for each np 
scattering angle bin by the fraction of high energy-loss 
tail events that survive the subtraction in AE pulse- 
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FIG. 17: The reconstructed (a) ytag spectrum for the CH2 tar- 
get before background subtraction, and (b) AE pulse-height 
spectrum for CH2-C for the 9p'^ bin 3-6°. The superimposed 
curves in both frames represent fits used in estimating back- 
ground subtraction accuracy. 



height spectra (see Fig.lSfc)). The tail events were inte- 
grated by summing all counts at AE values more than 
4cr above the center of a Gaussian fitted to the free- 
scattering peak, as illustrated in Fig. irTT b'). For the three 
largest 9^'^ bins studied, this approach breaks down be- 
cause the free-scattering AE peak develops a substan- 
tial Landau tail and is no longer well reproduced by a 
Gaussian shape. But for all (12) smaller-angle bins, the 
ratio of tail events after to before C subtraction fluctu- 
ates about zero, with a weighted average over angle bins 
of (2.97 ± 0.24) X 10~^. This measure is sensitive to the 
target thickness ratio as well as to the relative flux nor- 
malization for CH2/C. The two measures combined do 
not give compelling evidence of a need for any correc- 
tion, and are conservatively summarized by associating 
with the C subtraction an angle-independent correction 
factor C6 = 1.000 ± 0.004. 

The uncertainty estimated in this way also subsumes 
two other potential sources of systematic error. One is 
accidental coincidences between a real tagged neutron 
and an uncorrelated forward-going proton emerging from 
the GJT or the Cooler beam pipe (the most abundant 
sources of protons). To the extent that such coincidences 
passed all our cuts, they might have contributed to N2{0) 
for the CH2 target. However, since these accidentals are 
independent of the presence or nature of the secondary 
target, they would be subtracted via the C target mea- 
surements. The second effect concerns possible proton 
attenuation before the hodoscope, which is required as 
part of the event stream 2 hardware trigger. The dom- 
inant material between secondary target and hodoscope 
that might have served as a source of such proton losses 
is the AE scintillator, where tertiary interactions should 
cause abnormal energy loss. Since the cg uncertainty esti- 
mate includes allowance for such abnormal pulse heights 
in carbon-subtracted AE spectra, it should also include 
such proton attenuation effects. 

One further potential complication with the back- 
ground subtraction could have arisen if there had been 
any appreciable hydrogen buildup on the graphite tar- 



get used for the subtraction, a possibility Hmited by the 
hydrophobic nature of graphite. In this circumstance 
we would subtract some small fraction of the valid free- 
scattering events, and thus would introduce an effective 
overall normalization error in the hydrogen target thick- 
ness tn used in Eq. ^ To estimate this effect we consid- 
ered np scattering events forward of Oc.m. = 90°, where 
event stream 2 contained some coincidence events, with 
a forward-scattered neutron detected in the hodoscope 
and the larger-angle proton detected in the AE scintil- 
lator and (at least) the front two MWPCs. The angle 
of the proton was determined from MWPC ray-tracing, 
while that of the neutron was deduced from the ho- 
doscope elements fired and from the position inferred 
from the time difference between hodoscope phototubes 
mounted at the two ends of each element 0. The 
opening angle spectrum reconstructed for such np coin- 
cidence events exhibited a clear free-scattering kinematic 
peak for the CH2 target, but only the Fermi-smeared 
and acceptance-limited angular correlation characteris- 
tic of quasifree scattering for the C target (see Fig. I18|l . 
Figure ITSl includes fits to the distributions for both tar- 
gets based on the sum of a quadratic background and 
a Gaussian free-scattering peak, with the peak location 
and width fixed for the C target to the values determined 
from CH2. The fit for C is statistically consistent with no 
hydrogen content in the graphite target, with a Icr limit 
on the hydrogen thickness of 0.4% that of the hydrogen 
in CII2. We thus apply a cross section correction factor 
C7 = 1.000 ± 0.004 for hydrogen in the C target. 



2. Uncertainty in subtraction of corrupted events 

For the 2-stop event stream, we followed the procedure 
described in Sec. IIII Bl to subtract the punch-through 
events that had been corrupted by the electronic loss 
of backing detector pulse height information. There is 
no evidence for any systematic deviation in distribution 
shapes between the corrupted sample and our simula- 
tion of this sample using valid recorded punch-through 
events. Thus, the only uncertainty we consider is that 
in the normalization of the simulated sample to the cor- 
rupted events in stream 2. The normalization factors 
were determined from fits for the subsample of corrupted 
events that had valid backing detector timing signals, and 
the uncertainty in these normalization factors was then 
deduced from the change in normalization factor that 
caused an increase of unity in the overall value for the 
fit. The effect of this normalization uncertainty on the 
extracted 2-stop cross sections was typically « 0.01%, 
and is negligible in comparison with other systematic er- 
rors. Hence, we assign a correction factor cs — 1.000 with 
uncertainty < ±0.001 to the subtraction of corrupted 
events. 



16 



= 400 



100 







1 np Coincidences 


■ 




'\ CH2 Target 





















60 70 80 90 100 110 120 





1 ' 


1 ' 1 ' 1 ' 1 ' 

np Coincidences 
C Target 










1 . 


♦Hi 

1.1. . . . . 



Ol , 1 , 1 , 1 , l>^ .... j 

60 70 80 90 100 110 120 

Opening Angle (deg.) 



FIG. 18: np scattering opening angle spectra reconstructed 
for events where a large-angle proton fires at least the first 
two MWPCs and a forward neutron appears to fire the rear 
scintillator hodoscope. For such coincidence events, a clear 
free-scattering peak is seen with the CH2 production target 
(upper frame), while no hint of such a peak is seen for the C 
target (lower). The solid curves are fits with a Gaussian peak 
superimposed on a quadratic background. The distribution 
shape for C reflects the quasifree np scattering opening-angle 
spectrum convoluted with the coincidence acceptance of the 
forward detector array. 



3. Losses via software cuts 

The efficiency of software cuts applied to event stream 
2, but not to streams 1 and 3, was judged by comparing 
the ratio of cross sections obtained, after CH2 - C sub- 
traction, for all events failing vs. satisfying a given cut. 
The most important of these cuts were on AE{9p'^) (see 
Fig. and on Xtrack - Xtag and y track - Vtag - The lat- 
ter cuts were already dealt with in Sec. IV A 51 (Another 
cut, on the quality of track fits, is treated together with 
wire chamber inefficiencies below.) The Ai? cut limits 
were somewhat tighter than the 4<t allowance used in esti- 
mating background subtraction accuracy (see Sec. lVB T|l . 
We found the ratio of background-subtracted events fail- 
ing/satisfying the /S.E cut to be 1.0%, averaged over 
all event streams and angles. There is no evidence for 
any significant angle-dependence in this loss, but there 
are strong enough fluctuations in the losses from an- 
gle to angle or event stream to event stream that we 
assign a ±50% uncertainty to the losses. We thus ap- 



ply a corresponding, angle-independent correction factor 
cg — 1.010±0.005. With this systematic uncertainty, ap- 
plication of the AS cut still reduced the overall cross sec- 
tion error bars slightly because the quasifree background 
to be subtracted decreased significantly. 

4-. Reaction tail losses beneath the hodoscope energy 
threshold 

Since we did not use any software cuts on energy de- 
position in the rear hodoscope, we avoided the large 
corrections that would have been needed to account for 
protons lost to nuclear reactions in this hodoscope (see 
Fig. I10|l . However, if the reaction is sufficiently severe 
that the deposited energy falls below the hodoscope hard- 
ware threshold, then the event will have been lost in hard- 
ware to a trigger inefficiency. To estimate these potential 
losses, we fit the hodoscope energy spectra after CH2 - C 
subtraction to the sum of a Gaussian and an exponential 
(reaction) tail, as shown in Fig. EH The tail was extrapo- 
lated to zero energy deposition, and the ratio of yields be- 
low to above threshold (typically set at 5-10 MeV proton 
energy) was thereby estimated. The loss below thresh- 
old was found to be quite consistent with 0.4% for each 
scattering angle bin, so that we again have applied an 
angle-independent correction factor cio — 1.004 ± 0.002. 

5. Neutron polarization effects 

While the stored proton beam in the Cooler was un- 
polarized for this experiment, the neutron production 
reaction selected neutrons scattered to one side of the 
beam (beam right) at about 14° in the laboratory frame. 
At this angle, the D(p,n) charge exchange reaction that 
dominates our tagged beam production has a small po- 
larization, so that the beam neutrons would have been 
slightly polarized vertically (perpendicular to the hori- 
zontal production plane). The magnitude of this effect 
is PP'""'' w —0.1, where the minus sign indicates that 
for neutron production to the right of the cooled proton 
beam, the neutron spin points preferentially downward 
at the secondary target. The tagged neutron polariza- 
tion can then give rise to a left-right asymmetry in np 
scattering events: 

e^p{eA) = Pr''Anp{e)cos{^). (2) 

Measurements and phase shift solutions at intermediate 
energies show the np scattering analyzing power, 
A„p, to be a strong function of scattering angle, but with 
magnitude no larger than 0.12 over the angle range of 
interest for the present experiment. The cos(0) factor 
reflects the fact that it is only the component of the ver- 
tical neutron polarization perpendicular to the scattering 
plane for a given np event that matters. 

Since the scattering yield is simply redistributed be- 
tween scattering toward the left and the right, there 
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would be no effect at all on cross sections measured with 
a fully left-right symmetric forward detector array. Thus, 
the only residual polarization effect changes the measured 
yield by a fraction: 

S{e)=Pr''Ar,p{9) r a{e,(p)cos{4>)d4>l r a{e,cj,)dcf>, 

JQ Jo 

(3) 

where a{9, (p) is the fractional detector acceptance (deter- 
mined from fits such as those in Fig. Ill|l in the specified 
angle bin. The sign convention used here is that positive 
d{9) implies that we observe a higher event stream 2 yield 
than we should in the corresponding 9^^^ bin, necessitat- 
ing a correction factor cii{9p'^) — 1.0 — S{9). 
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FIG. 19: The estimated fractional cross section error intro- 
duced by neutron polarization effects for the 1-punch and 2- 
stop event samples, plotted as a function of np scattering 
angle. 

Figure [H^ shows the S{9) distribution calculated for 1- 
punch and 2-stop data samples from Eq.|3 taking Anp{9) 
from Nijmegen PWA93 calculations 24]. We find 6 = 
for both samples at all angles 9p'^ < 25°, because the de- 
tector has full azimuthal acceptance in that region. The 
small corrections have opposite sign, and hence tend to 
cancel, for the two samples at larger angles because the 2- 
stop events preferentially populate the left side of the sec- 
ondary target, while the 1-punch events originate mostly 
on the right (see Fig.Q. The latter difference is reflected 
in their respective a{9,(f)) functions (see Fig. Ill|l used in 
Eq.El 

While the np analyzing power and the forward detector 
acceptance functions are well determined in this experi- 
ment, we assign a significant uncertainty to the average 
tagged neutron polarization, = -0.10 ± 0.05, to 

account for contributions from production mechanisms 
other than charge exchange. There is correspondingly 
a ±50% uncertainty assigned to each value of 6{9) in 
Fig. El but these errors are completely correlated from 
one angle bin to another and are strongly correlated be- 
tween 2-stop and 1-punch data samples. The largest net 



uncertainty from neutron polarization after the (sepa- 
rately corrected) 1-punch and 2-stop results are combined 
is ±0.6 X 10~^, and so we conservatively assign an angle- 
independent uncertainty of ±0.001 to cn. 



C. Target thickness, acceptance and efficiency 
errors 

1. CH2 target thickness uncertainty 

The overall normalization uncertainty associated with 
tn in Eq. n was determined to be ±0.4% (i.e., C12 = 
1.000 ± 0.004) from careful weighing of the CH2 target 
used. Because the target only sat in a secondary neutron 
beam of low flux, there should not have been any appre- 
ciable deterioration in the target during the length of the 
run, nor was any visually evident. 



2. Acceptance uncertainty 

The acceptance uncertainty was determined indepen- 
dently for each angle bin by varying the most critical one 
or two detector geometry parameters used in the simu- 
lations (see Fig. Illf) from their best-fit values until the 
overall value for the simulated vs. measured (jj'^p dis- 
tribution in that angle bin increased by unity. Because 
the optimized values of per degree of freedom for the 
different angle bins and event samples were statistically 
distributed about 1.14, rather than 1.00, we multiplied 
these acceptance changes by a uniform factor of 1.07 
to arrive at final systematic uncertainties. The accep- 
tance uncertainty is strongly angle-dependent, varying 
from ±0.001 at 9c.m. > 120°, where > 95%, to ±0.017 
at 9c. m. = 90°, where « 50%. With this evaluation 
method, we consider the estimated uncertainties to be 
largely uncorrelated from angle bin to bin. 



3. Wire chamber efficiencies 

The MWPC's were not used at all in forming a hard- 
ware trigger, but in software event reconstruction we re- 
quired at least one hit registered in each of the x and 
y planes, for chambers 1, 2 and 3, plus a value be- 
low an upper threshold for fitting these hits to a straight 
line track. Thus, the overall MWPC efficiency to use in 
extracting absolute cross sections is: 

miwpc = Vxil)Vyi'^)Vx{'2)rjy{2)ri^{3)rjy{3)ri^^ quality 

(4) 

The efficiency of each MWPC plane was determined from 
tracks reconstructed without the benefit of the plane in 
question, based on the fraction of such tracks that pro- 
duced a hit on this plane in the immediate vicinity of 
the reconstructed crossing point. Each of the first six 
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factors in Eq.0|was found to exceed 0.99 and was deter- 
mined with an uncertainty « ±5 x 10~^. Their product 
is 0.985 ±0.0013. 

The factor Tyg^ quality ~ 0.998 ± 0.0011 was deter- 
mined by estimating the number of free-scattering events 
removed from the analysis by the cut. This was done 
by examining IS.E spectra for individual angle bins, fol- 
lowing carbon subtraction, for events that failed the fit 
quality test. The number of free-scattering events (and 
its uncertainty) in each such spectrum was extracted by 
fitting Gaussian peaks of the same position and width as 
those used for the normal AE spectra, such as Fig. llTf b). 
There was no indication in these analyses that any of the 
factors in Eq.^] varied with position on the MWPC's, or 
therefore with scattering angle. The overall wire cham- 
ber efficiency correction is thus an angle-independent 
ci4 = 1.0017 ±0.002. 



4. Trigger inefficiencies 

Inefficiencies in detectors used to form the hardware 
trigger lead to loss of events in an unrecoverable way. 
Possible tagger inefficiencies do not matter here, because 
they lead to loss of the same fraction of events from 
streams 1, 2 and 3, and hence do not affect the cross 
sections determined from ratios of event yields in these 
streams. The two detectors used to form the hardware 
trigger for event stream 2, but not for stream 1, are the 
AE scintillator and the rear hodoscope. The former was 
viewed by four phototubes, at least three of which were 
required to give signals surpassing threshold in the trigger 
logic. In the data analysis, we were able to determine for 
each scattering angle bin the ratio of reconstructed free- 
scattering events that had only three vs. all four AE pho- 
totubes above threshold. We then estimated the /S.E trig- 
ger inefficiency under the conservative assumption that 
the ratio of free-scattering events with two or fewer pho- 
totubes firing to those with three firing would be the 
same as the determined ratio of events with three to four 
firing. (Some illuminated locations on the scintillator 
lacked a direct line of sight to one or another, but never 
simultaneously to two, of the four phototubes.) The re- 
sulting inefficiency appears to show a systematic angle- 
dependence, roughly represented by 0.008 cos^(^p''), i.e., 
the inefficiency grows as the AE pulse height shrinks. 

We have considered two different types of potential 
hodoscope trigger inefficiencies. Problems in an individ- 
ual hodoscope element or phototube would show up as 
an inefficiency localized in 9 and (j), and therefore as a 
deviation of the measured ^-distribution for some an- 
gle bins from the simulated acceptance function. Any 
such localized trigger-level inefficiencies should thus be 
subsumed in the acceptance uncertainty calculation men- 
tioned above. 

However, an electronic inefficiency in the modules 
forming the hodoscope trigger logic could have caused 
equal fractional losses in all angle bins. A limit on this 



inefficiency was estimated from event stream 4 (observ- 
ing tagged protons from the GJT), which included the 
Veto2 scintillator directly in front of the hodoscope, but 
not the hodoscope itself, in the trigger logic. We found 
that (0.6 ± 0.1)% of these triggered events were not ac- 
companied by hodoscope signals above threshold in both 
relevant phototubes, of which 0.4% have already been 
accounted for as reaction tail losses below threshold (see 
Sec.|VB4j. 

Combining the above effects, the overall correction fac- 
tor for trigger inefficiencies has been taken as C15 = 
[1.002+0.008 X cos2(6i^=)] ± [0.001+0.004 x cos^iO^j,")] . The 
angle-dependent part of the uncertainty here is intended 
to accommodate observed fluctuations in the inferred A_E 
trigger efficiency, and is viewed as largely uncorrelated 
among different angle bins. 

5. Dead time differences among event streams 

Triggers in all event streams were blocked electroni- 
cally at an early stage in the event trigger logic by a 
common busy signal reflecting electronic readout or com- 
puter processing activity in any of the event streams. To 
first order, then, the different streams should have a com- 
mon dead time (« 10% for typical running conditions), 
and the dead time should cancel in the event stream ra- 
tios from which cross sections are deduced (see Eq. . 
However, this cancellation is imperfect, as revealed by ra- 
tios of scaler values recording the number of tagged neu- 
tron vs. tagged proton candidates before and after busy- 
vetoing. Typically, ~ 1% fewer neutron tags survived the 
veto, and this was traced to the occurrence of bursts of 
electronic noise triggers from the tagger. While the loss 
of these noise triggers should not have directly depleted 
the valid sample of any event streams (i.e., (2)-(4)) that 
required other detectors in coincidence with the tagger, 
it did reduce the number of valid events recorded in the 
neutron flux stream (1), because all raw neutron tags, 
whether valid or not, contributed equally to the count- 
down of a (divide by 20) prescaler used for this stream. 
To compensate for this loss of neutron flux events, we 
have to reduce the extracted cross sections at all angles 
by a factor cig = 0.991 ± 0.005. The uncertainty in this 
correction allows for possible model-dependence in our 
interpretation of the live-time difference inferred from the 
scaler ratios. 



D. Errors in Determination of Kinematic Variables 
for np Scattering 

1. Neutron energy errors 

As explained in Sec. lIVI the data were analyzed in nar- 
row neutron energy slices, with each result then being 
corrected slightly in order to extract a final overall cross 
section at the single mean energy of 194.0 MeV. There is 
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an overall scale uncertainty in the tagged neutron ener- 
gies that we estimate to be ±150 kcV, with roughly equal 
contributions from the energy of the stored primary pro- 
ton beam in the Cooler and the energies extracted from 
the tagger for the low-energy recoil protons. The stored 
beam energy (202.46 MeV) is based on the precisely mea- 
sured rf frequency (1.96502 MHz) and the Cooler circum- 
ference, which has been previously calibrated to bet- 
ter than 1 cm out of 87 m, translating to ±70 keV. In 
the "coasting" (rf off) mode used for data taking, the 
beam energy is maintained by interactions with the cool- 
ing electrons, and this may increase the beam energy un- 
certainty to ~ 100 keV. The energy scale of the recoil pro- 
tons is calibrated by analysis of ^^^Th a-source spectra 
measured with the tagger ^ind its ±100 keV uncer- 
tainty arises predominantly from thickness uncertainties 
for detector dead layers, combined with the quite differ- 
ent corrections for energy loss in the dead layers needed 
for protons vs. the calibration a-particles. 

The energy scale uncertainty could be translated into 
a consequent cross section uncertainty as a function of 
angle by using Nijmegen PWA93 calculations to evaluate 

So energy {01") = (±150 kcV) — MoV- (5) 

Although this systematic error can be angle- 
dependent, the values at different angles would still be 
completely correlated, since the neutron energy scale will 
be off in the same direction for all angles. Hence, we pre- 
fer to not include this effect in the cross section uncer- 
tainties, but rather to quote the measured cross sections 
as applying at a mean neutron energy of 194.0 ± 0.15 
MeV. 



tag, y tag) of the Scattering origin on the secondary tar- 
get, and over all scattering angles. Consistent values were 
extracted for the 2-stop and 1-punch samples. 

The angle uncertainties were estimated within 1x1 
cm^ pixels in {xtag, Utag) as half the mean event-by-event 
difference between angles reconstructed by two different 
approaches. In the case of Oinc one method utilized tagger 
information only to predict the neutron trajectory, while 
the second considered instead the straight line from the 
neutron production vertex on the GJT, inferred from the 
tagger, to the intersection (xtrackiVtrack) of the recon- 
structed forward proton track with the secondary tar- 
get. For 9 out we used proton tracks reconstructed with 
MWPC geometry parameters that were either (i) op- 
timized to minimize the overall value for tracks, or 
(ii) adjusted to increase overall by unity. A yield- 
weighted average of the results over all target pixels gives 
{Sdinc) = 1-3 mrad and (SOout) = 0.04 mrad. 

The cross sections were not corrected for potential sys- 
tematic angle errors, but we extract from Eq. Elnet sys- 
tematic uncertainties of ±0.4% for 120 < 6c.m. < 180°, 
±0.3% for 100 < dc.m. < 120° and ±0.1% for 90 < 
0c.m. < 100°. Because the extracted incidence angle 
differences (between the two methods described above) 
exhibit sizable fluctuations from one target pixel to an- 
other, or from one angle bin to another, we view these 
estimated uncertainties as uncorrelated from angle bin to 
angle bin. 



E. Summary of angle-dependence 



2. Scattering angle errors 

A systematic error 56'^^ in determination of the cen- 
troid np scattering angle within a given analyzed bin is 
equivalent to an error Saangie{Sp'^) in the measured dif- 
ferential cross section: 

^^d[da/dn]{0'p'') 

Scrangleidp"^) = SOp'^ |l94MoV, (6) 

where the angular derivative of the cross section can be 
taken, for example, from Nijmegen PWA calculations 
[23 |. In evaluating SOp"^, we consider the contributions 
from uncertainties SOinc in the neutron incidence angle 
on target deduced from the tagger and 60out in the angle 
of the outgoing proton through the MWPC's: 



se;'^ = m^c)' + {S9out)']'^', (7) 

where the averages are evaluated over the full free scat- 
tering event sample, over the transverse coordinates 



The effect of the correction factors Ci associated with 
the various sources of systematic error considered in this 
section is cumulative, and averages 1.10, with small vari- 
ations with angle and data sample, as summarized in 
Table We assume, however, that the various uncer- 
tainties are uncorrelated with one another, and we add 
them in quadrature to obtain final systematic error esti- 
mates. The majority of error sources we have considered 
are explicitly or effectively angle-independent; when com- 
bined, these yield an overall normalization uncertainty 
of ±1.5%. The uncertainties associated with our mea- 
surements of acceptance and scattering angle, and with 
trigger inefficiencies, are considered angle-dependent and 
uncorrelated from point to point. These three sources 
are combined to give the net point-to-point systematic 
uncertainties in Table Ih] where we also collect our fi- 
nal absolute cross section measurements obtained from 
a weighted average over the three independently ana- 
lyzed and corrected data samples (1-punch, 2-stop with 
Epi < 5.0 MeV and 2-stop with Epi > 5.0 MeV). The 
point-to-point and normalization uncertainties combine 
to give an overall systematic error of ±1.6% in most an- 
gle bins. 
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TABLE II: Final differential cross section results for np scat- 



tering at 


= 194.0±0.15 MeV, 


averaged over data samples. 


cm. angle 


{daldiL)c.m. 


Stat. unc. 


Syst. unc." 


(deg.) 


(mb/sr) 


(mb/sr) 


(mb/sr) 


92.7 


1.98 


0.06 


0.03 


98.8 


2.00 


0.05 


0.02 


104.8 


2.31 


0.05 


0.02 


110.8 


2.57 


0.05 


0.02 


116.8 


3.01 


0.05 


0.02 


122.8 


3.47 


0.06 


0.02 


128.8 


4.01 


0.06 


0.02 


134.9 


4.75 


0.07 


0.03 


140.9 


5.34 


0.08 


0.03 


146.9 


5.98 


0.08 


0.04 


152.9 


6.64 


0.10 


0.04 


159.0 


7.61 


0.11 


0.05 


165.0 


8.89 


0.14 


0.06 


171.0 


10.62 


0.19 


0.07 


177.0 


11.86 


0.34 


0.08 



"This column lists point-to-point systematic uncertainties. In ad- 
dition, there is an overall cross section scale uncertainty of ±1.5%. 



VI. DISCUSSION 

The preceding section provided a detailed catalogue 
of the issues that must be carefully controlled to mea- 
sure precise absolute cross sections with medium-energy 
neutron beams. To our knowledge, no previous exper- 
iments have attempted a comparable degree of control. 
The best existing absolute neutron-induced cross section 
standards at intermediate energies are from attenuation 
measurements of total cross sections [27], which are not 
suitable for calibrating neutron fluxes. It is hoped that 
the present results will provide a new calibration stan- 
dard. The excellent agreement of our experimentally 
determined absolute cross section scale with that given 
by the Nijmegen PWA93 solution (see Fig. I14|) confirms 
the consistency of our results with the total cross section 
measurements. 

The level of agreement of our measurements with 
PWA's at En — 194 MeV is presented in more detail 
in Fig.|2Dl Here it is seen that, while the absolute cross 
section scale of the experimental results is in very good 
agreement with the Nijmegen PWA93 solution, there is 
a small systematic deviation in angular shape between 
the two: our results are higher than PWA93 by 2-3% 
for 135 < Oc.m. < 165° and lower by a similar amount 
for 100 < Qc.m. < 130°. These deviations consider- 
ably exceed our estimated systematic uncertainty in the 
angle-dependence. In particular, we note that the for- 
ward detector acceptance used in the former angle range 
is already essentially 100% (see Fig. , so that the ex- 
tracted cross section cannot be overestimated by virtue 
of underestimating acceptance. Furthermore, the results 
for the three independently analyzed data samples agree 
extremely well in this angle region (see Fig. I13|l. We do 
see a possible small, statistically marginal, systematic de- 



viation among our three data samples in Fig.^|over the 
angle region from 100 to 130°, with the 1-punch cross 
sections falling on average a few % below those for the 
two 2-stop event samples. However, even if this difference 
reflects a real systematic problem, it could only pull the 
averaged cross section down by less than 0.5% in this re- 
gion, too small to account for the deviation from PWA93 
in Fig. 
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FIG. 20: The relative differences of the present absolute np 
scattering differential cross sections and of two SAID PWA 
solutions 28] from the Nijmegen PWA93 solution all 
at £;„ = 194 MeV. The SP40 solution is from a 2003 analysis 
of the database from 0-400 MeV, while SP05 is the current 
SAID solution, fitted over the range 0-3.0 GeV, including the 
present data in the fit. 



We also show in Fig. [201 the relative differential cross 
section differences between two recent SAID PWA so- 
lutions 0, Hi] and the Nijmegen PWA93. The various 
PWA solutions differ from one another by as much as 
2-3% also in the angle region displayed. Furthermore, 
we note that the SAID solution has shifted by ~2% after 
inclusion of the present results in the fitted np database 
(even though that inclusion was carried out by adding our 
full, mostly angle-independent, systematic uncertainties 
in quadrature with our statistical uncertainties, thereby 
under- weighting the present data in the fit). We con- 
clude that the deviations between the present results and 
PWA93 are of the same order as the present uncertainties 
in the PWA solutions, and most likely point to the need 
to refit phase shifts. Wc note, however, that there is a 
conceptual flaw in the procedures for such refitting to a 
database where all experiments have systematic uncer- 
tainties, but there is considerable variability in the level 
at which those systematic uncertainties are reported in 
the literature. 
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Finally, we address the comparison of the present re- 
sults with those from the recent experiments by the Upp- 
sala and Freiburg groups, both of which have been 
rejected from the np database used in the Nijmegen and 
SAID PWA's. As illustrated in Fig.^^by the comparison 
of the two experimental results with PWA curves at the 
respective energies of the experiments, the present results 
deviate systematically from those of in the steepness of 
the back-angle rise in cross section. These deviations are 
larger than the differences anticipated from the difference 
in neutron energy between the two experiments. There 
is a similar, though not quite as pronounced, systematic 
deviation of the present results from those of Franz, et 
al. 0, shown in Fig.Q] 

It is difficult to say definitively whether there might be 
a common problem that caused excessive cross sections 
near Oc.m. — 180° in both of these earlier, conipletely 
independent and quite different, experiments [3,I3- We 
note only that measurements near 9^'^ = 0°, where the 
solid angle is vanishing, can be tricky with a secondary 
neutron beam of sizable angular divergence. One of the 
great advantages of the use of a tagged beam is that 
we are able in the present experiment to determine the 
neutron incidence angle event by event. We thereby see 
that an angle bin near 9^^^ = 0° with respect to the 
central neutron direction in fact includes important con- 
tributions from 9p^ > (so that \d cos{9p'^)\ for the bin 
in Eq. ^exceeds |cicos(6'p°^)|), arising from neutrons de- 
viating from the central beam direction. In Fig. 1211 we 
compare the present results with those we would have ex- 
tracted if we had chosen to ignore the neutron incidence 
angle information from the tagger in reconstructing the 
np scattering angle event by event. Such neglect is seen to 
give rise to a systematic overestimate of the cross section 
at the largest angles by ^ 5%, comparable to the system- 
atic deviations of the results in Refs. 0, |^ from those of 
the present experiment. (It also leads to a substantial 
underestimate near 9c.m. — 90°, where the yield is falling 
rapidly with increasing 6'^"'', due to detector acceptance 
edges.) The effect would differ for different experiments, 
depending on the angular profile of the neutron beam, 
including any effects from scattering off collimator edges 
(the present experiment used no collimators). 



VII. CONCLUSIONS 

A tagged intermediate-energy neutron beam produced 
at the lUCF Cooler ring has facilitated a measurement 
of the np scattering differential cross section at 194 MeV 
bombarding energy to an absolute precision « ±1.5% 
over the cm. angular range 90°-180°. The usage of 
carefully matched and frequently interchanged solid CII2 
and C secondary targets permitted an accurate back- 
ground subtraction, reducing reliance on kinematic cuts 
that might have introduced larger systematic uncertain- 
ties. The internal consistency in both magnitude and 
angular shape of the cross sections extracted from in- 
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FIG. 21: The effect on the present analysis of neglecting 
tagger information about the neutron incidence angle in the 
reconstruction of the np scattering angle for each event. 
The closed circles represent the final results, while the open 
squares are those when the scattering angle is estimated only 
with respect to the central neutron beam direction. Note the 
suppressed zero on the cross section scale. 



dependent data samples characterized by substantially 
different neutron beam spatial and energy profiles sup- 
ports the accuracy of the tagging technique. System- 
atic uncertainties in the measurement, affecting both the 
overall absolute scale of the cross sections and the angu- 
lar dependence, have been carefully delineated, often via 
auxiliary measurements and analyses. 

The present results are in reasonable agreement with 
the Nijmegen PWA93 calculation, over the full angular 
range covered, although there are systematic deviations 
at the 2-3% level in the angular dependence that might 
be removed by minor tuning of phase shifts. In contrast, 
the present results deviate systematically from other re- 
cent measurements , especially in the steepness of 
the back-angle cross section rise. Our results thus appear 
to validate the omission of these earlier experiments from 
the database used in partial-wave analyses of np elastic 
scattering, while also suggesting a possible experimen- 
tal cause of the earlier overestimates of the cross section 
near 9c.m. — 180°. As the back-angle rise is particularly 
influential in pole extrapolations that have occasionally 
been used [12. to extract the charged pion-nucleon- 
nucleon coupling constant the present data also bear 
on that coupling strength. Since our measurements at 
the largest angles are consistent with, or even slightly less 
steep than, the PWA93 solution, a valid pole extrapola- 
tion analysis of the present results should yield a coupling 
constant value no larger than that (/^ = 0.0748±0.0003) 
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extracted from the Nijmegen PWA 
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